Draft version February 2, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. ^jllj^A, 



UV LIGHT FROM YOUNG STARS IN GEMS QUASAR HOST GALAXIES AT 1.8 < Z < 2.75 



K. Jahnke . S. F. Sanchez , L. Wisotzki ■ , M. Barden"*, S. V. W. Beckwith •', E. F. Bell', A. Borch'' 



J. A. R. Caldwell^, B. Haussler"' 



C. Y. PENG^ H.-W. Rix-'* 



o 
o 

(N 
^3 



> 
(N 

m 
o 

o 
o 



S. JoGEE^, D. H. McIntosh^, K. Meisenheimer'' 
R. S. Somerville^ and C. Wolf^ 

Draft version February 2, 2008 

ABSTRACT 

We have performed ffST" imaging of a sample of 23 high-redshift (1.8 < z < 2.75) Active Galactic 
Nuclei, drawn from the COMBO- 17 survey. The sample contains moderately luminous quasars {Mb ~ 
—23). The data are part of the GEMS imaging survey that provides high resolution optical images 
obtained with the Advanced Camera for Surveys in two bands (F606W and F850LP), sampling the 
rest-frame UV flux of the targets. To deblend the AGN images into nuclear and resolved (host galaxy) 
components we use a PSF subtraction technique that is strictly conservative with respect to the flux 
of the host galaxy. We resolve the host galaxies in both filter bands in 9 of the 23 AGN, whereas 
the remaining 14 objects are considered non-detections, with upper limits of less than 5 % of the 
nuclear flux. However, when we coadd the unresolved AGN images into a single high signal-to-noise 
composite image we find again an unambiguously resolved host galaxy. The recovered host galaxies 
have apparent magnitudes of 23.0 < F606W < 26.0 and 22.5 < F850LP < 24.5 with rest-frame 
UV colours in the range —0.2 < (F606W — F850LP)obs < 2.3. The rest-frame absolute magnitudes 
at 200 nm are —20.0 < M200 nm < —22.2. The photometric properties of the composite host are 
consistent with the individual resolved host galaxies. We find that the UV colors of all host galaxies 
are substantially bluer than expected from an old population of stars with formation redshift z < 5, 
independent of the assumed metallicities. These UV colours and luminosities range up to the values 
found for Lyman-break galaxies (LBGs) at 2; = 3. Our results suggest either a recent starburst, of 
e.g. a few per cent of the total stellar mass and 100 Myrs before observation, with mass-fraction and 
age strongly degenerate, or the possibility that the detected UV emission may be due to young stars 
forming continuously. For the latter case we estimate star formation rates of typically ~6Mq yr~^ 
(uncorrected for internal dust attenuation), which again lies in the range of rates implied from the 
UV flux of LBGs. Our results agree with the recent discovery of enhanced blue stellar light in AGN 
hosts at lower redshifts. 

Subject headings: galaxies: active - galaxies: high-redshift - galaxies: fundamental parameters - 
galaxies: starburst - quasars: general 



1. INTRODUCTION 

Around redshifts of z ~ 2-3, luminous quasars were or- 
ders of magnitude more numerous than today. Although 
the physics of how active galactic nuclei evolve is still 
not understood, several links between galaxy and quasar 
evolution have emerged over recent years. The obser- 
vational confirmation of supermassive black holes in the 
nuclei of all gal axies with a substantial bulge component 
(e.g. [Gcbhardt et alJl2000H makes every such galaxy a 
potential AGN host. The strong evolution of the AGN 
space density could therefore be related to the availability 
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of accretion fuel in the host galaxies, or to the frequency 
of AGN triggering events. 

Gravitational interaction and major or minor merging 
of galaxies have long been suggested as important fac- 
tors in driving nuclear activity in galaxies. Confirming 
any of these as the dominant process has proved difficult, 
mainly because the morphological characteristics found 
for relatively nearby AGN host galaxies are so diverse. 
Furthermore, the properties of the hosts in the 'heyday' 
of quasars (z > 2) are still elusive, a consequence of 
the substantial observational difficulties. The contrast 
between the bright nuclear point sources and the sur- 
rounding galaxy increases dramatically beyond z 1 as 
a result from both surface brightness dimming and wave- 
band shifts towards the rest frame UV. 

The last years have seen numerous attempts to re- 
solve the host galaxies of high-redshift quasars. Ow- 
ing to the observational challenges of detecting dis- 
tant host galaxies the observational effort for each 
object is large, and the observed samples have con- 
sequently been very small, of the order of < 5 
per target group. While radio-loud quasars appear 
to be very e xtended and have been resolved out to 
z ~ 4 (e.g., ILe hnert et al." 199?; 'Carballo et al.' 1991 
Hutchings et al. 1999: Kukula ct al._2001; Hutching! 
2003HSancher fc Gonzalez-Serrancil2003|) . this is not the 
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case for the large majority of radio-quiet quasars. 

At high redshifts two constraints dominate observa- 
tional studies of host galaxies: On one side, very good 
seeing conditions are required to maximize the spatial 
contrast of the compact nuclear source compared to the 
extended host galaxy. On the other, large telescope 
apertures are preferrential to trace faint quasar hosts 
to as far away from the nucleus as possible. Thus sig- 
nificant progress had to wait for 8m-class telescopes at 
very good sites with active optics systems - with a very 
high light collecting power but atmospheric seeing limi- 
tations - and for the HST and its high space-based sen- 
sitivity, combined with unprecedented spatial resolution, 
but limited size that might miss light from faint outer 
structures of the hosts. Some host galaxies of radio- 
quiet quasars at 2; ~ 2 h ave now been resolve d both 
in the near infrared (Arctx aea et al.li l998: Kukul a et al.l 
l200lURidgwav et al.i2001i.Falomo et aL.20041 and in the 
optical domains ^Hutchings et al.l 120021) . showing these 
objects to be moderately luminous, corresponding to 
present-day L* or slightly brighter. 

However, host galaxy colours have been unavailable, 
precluding estimates of the mass-to- light ratio (M/L). 
Thus, without colours the observed luminosities, and 
their evolution with redshift, cannot be mapped to the 
mass evolution if young stars contribute a major frac- 
tion of AGN host's light. This is important as several 
high-luminosity quasars at z > 2 a ppear to be located 
in very UV-luin inous host galaxies llLehnert et alJll992l: 
lAretxaga et"al] 119981 iHutchings et al.ll2002D . Also, at 
low redshifts there is a link between nuclear activity 
and enhanced global star f ormation in the host galax- 
ies. iKanffmann et alJ 1)200,'^ reported that SDSS spectra 
of local Seyfert 2 galaxies show a significant contribution 
from young stellar populations, and that this trend is 
strongly correlated with nuclear luminosity. In a multi- 
colour study of low-z QSO hosts (Jahnke et al. 2004) as 
well as at intermediate redshifts (,Sanchez et al...2004. see 
below) we found that hosts of elliptical morphology can 
be significantly bluer than the bulk of inactive ellipticals. 
These results indicate that in the recent past the star 
formation activity in galaxies hosting an AGN may be 
different from normal galaxies. The details are far from 
understood. Clearly more information is required to in- 
vestigate the relation of starformation and AGN activity, 
their common cause or causal order and the evolution of 
these properties with redshift. 

The new generation of wide-field imaging mosaics ob- 
tained with the Hubble Space Telescope (HST), espe- 
cially in the conjunction with deep AGN surveys, has 
opened a new observational avenue towards AGN host 
galaxy studies. Here we present first res ults on AGN 
within the gems project ijRix et af]l2004j) . the largest 
HST colour mosaic to date. In the present paper we in- 
vestigate the presence of rest-frame ultraviolet light in 
a substantial sample of z > 1.8 AGN, all with nuclear 
luminosities near M b = — 23. In a companion paper 
ijSanche z et alJ I2004j) we study rest-frame colours and 
morphological properties of a sample of intermediate- 
redshift {z < 1) AGN. 

The paper is organised as follows. We first describe 
the sample selection and properties together with a sum- 
mary of the observational data (Sect.|2). We then com- 
ment on the decomposition of the nuclear and galaxy 



contribution, including a brief summary of the exten- 
sive simulations that we use to estimate measurement 
errors (Sect. In Sect. ^ we present the measured 
host galaxy magnitudes and describe our treatment of 
non-detections. We move on to discuss the results in 
Sect. 13 followed by our conclusions in Sect. El We 
use Ho = 70kms~iMpc~\ fi^ = 0.3 and Ha = 0.7 
throughout this paper. All quoted magnitudes are zero- 
pointed to the AB system with ZPpgoew = 26.493 and 
ZPF850LP - 24.843. 

2. AGN IN THE GEMS SURVEY 
2.1. Overall survey properties 

GEMS, G alaxy Evolution from Morphologies and SEDs 
l)Rix et alJ i2004') is a large imaging survey in two bands 
(F606W and F850LP) with the Advanced Camera for 
Surveys (ACS) aboard HST. Centered on the Chandra 
Deep Field South (CDFS), it covers an area of ~ 28' x 28' 
(78 ACS fields). Each ACS field was integrated for 
3 X 12-13 min exposures per filter (one orbit), dithered 
by 3" between exposures. The individual images were 
then combined, corrected for the ACS geometric distor- 
tion, and at the same time rebinned to a finer pixel grid 
of 0'.'03, achieving approximate Nyquist-sampling of the 
PSF. The image combination also removed artefacts such 
as cosmic ray hits and satellite trails. Resulting point 
source limiting magnitudes are mAB(F606W) = 28.3 
(5cr) and mAB(F850LP) = 27.1 (5cr). In its central ~l/5, 
GEMS incorporates the epoch 1 data from the GOODS 
ijGiavalisco et al.ll2004|) project, that are similarly deep 
as the other GEMS fields. Further details of the data re- 
duction procedure will be given in a forthcoming paper 
(CaldweU et al. 2004, in prep.). 

The area covered by gems coincides w ith one of 
four fields covered by the combo- 17 survey ijWolf et alJ 
|2004|) . which produced a low- resolution spectrophoto- 
metric data base (based on photometry in 17 filters) for 
about 10 000 galaxies and 60 type 1 AGN br ighter than 
R < 24 (Vega zeropoi nt) in the CDFS area l|Wolf et alJ 
l2003t iWolf 613111200^ . The large number of filters per- 
mitted simultaneous assignment of accurate SEDs and 
redshifts for both galaxies and type 1 AGN. Galaxies and 
AGN are classified by matching an SED template library 
to the set of 17 photometric points. The AGN SED is 
composed of a range of continuum spectra with added 
broad emission lines (all details are given in IWolf et alJ 
I2004|) . Type 2 AGN as well as very low luminosity AGN 
are invariably classified as galaxies. COMBO-17 photo- 
metric redshifts are very reliable, with an rms scatter 
of 0-2/(1 -f z) ~ 0.02 for galaxies (at z < 1.2) and 
Cz/ {1 + z) ~ 0.03 for AGN at all redshifts. In this paper 
we address specific COMBO- 17 sources just by their run- 
ning identifiers; the full COMBO-17 list of classifications 
in the CDFS will be made available in the future (Wolf 
et al. 2004, in prep.). 

2.2. The AGN sample 

The COMBO- 17-selected AGN in the gems field range 
over redshifts from z ~ 0.5 up to z ~ 4. In this study 
we investigate the high-z part of the AGN distribution. 
Our sample contains all AGN brighter than i? = 24 in 
the redshift range 1.8 < z < 2.75 that show a meaningful 
counterpart in the gems images. This excludes three 
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2.2 2.4 
redshift z 

Fig. 1. — COMBO-17 redshifts and R magnitudes of the AGN 
in the sample. The dashed line marks the expected R band mag- 
nitude for an AGN of Mg = — 23 at the corresp onding redshifts, 
assuming a typical quasar spectrum tWolfet aLlEOOS) , the dotted 
lines correspond to Mg = —21.5 and Mg = —24.5 for reference. 



3.1. Background and variances 

Even though space based, ACS shows a non-neghgible 
background from stray hght. In the reduction process 
already a global, outlier clipped median background was 
subtracted (Caldwell et al. 2004, in prep.). As the de- 
blending of nuclear and host galaxy component with two- 
dimensional modelling is sensitive to background sources, 
we applied an extra procedure to remove net residuals in 
the local background. This included an iterative mask- 
ing of all objects in the field and the determination of 
the local background from the object-free regions. For 
each square of 200 x 200 pixels an average from the un- 
masked pixels was computed, with a subsequent bilin- 
ear interpolation between these values to yield a back- 
ground estimate for the whole field. After background 
adjustment, small subimages of 128 x 128 pixels were ex- 
tracted around all AGNs, corresponding to a field of view 
of 3'.'84 X 3'.'84. This field size contains > 99 % of the ACS 
PSF flux. Since the AGN are not strongly resolved this 
fraction also applies to the AGN flux. 

The data reduction procedure kept record of the in- 
dividual pixel weights throughout the process of reduc- 
tion and combination. This information was then used in 
combination with the shot noise derived from pixel count 
rates to construct variance images which were later used 
in the error budget calculations. 



objects with a primary classification as AGN by COMBO- 
17 that were apparently low redshift (z ~ 0.1) emission 
line galaxies. 

At these redshifts, all light detected in the gems bands 
will originate in the rest frame ultraviolet. In fact, the 
lower redshift limit has been imposed to ensure that the 
long wavelength cutoff of the F850LP filter is still lo- 
cated below the Balmer jump for all objects. The up- 
per redshift boundary, on the other hand, was set to 
avoid contamination from possible extended Lya emis- 
sion in the F606W filter. The resulting sample contains 
23 AGNs, three of which are positioned in the overlap 
region of two tiles so two separate images exist for these. 
One object (combo-17 05696) shows inconsistent pho- 
tometry between the two tiles due to variability over the 
111 days between the two integrations. Another object 
(combo- 17 19965) was classified in combo-17 with a 
redshift of z = 0.634 that had to be revised following 
spectroscopy observations (G. Worseck, private commu- 
nication) and now entered this sample with z — 1.90. 
Table Q gives an overview over the sample properties, 
and Figure n shows the distribution in R and z. The av- 
erage absolute magnitudes for these objects - which we 
call intermediately luminous quasars - place them close 
to the canonical division of Mb — —23 between Seyfert 
galaxies (low-luminosity) and QSOs (high-luminosity). 

Several X-ray sources in the CDFS hav e already been 
studied by HST with the WFPC2 camera CSchreie r et"an 
[2001; Koekcmocret al. 2002; Grogin et al. 2003), most 
of them faint AGN, but their sample is completely dis- 
joint from the COMBO-17 AGN selection for this redshift 
range. One object falls into our redshift range but has 
R > 24, outside our selection limits. 

3. DATA ANALYSIS 



3.2. PSF estimation 

Compared to ground based telescopes, the point spread 
function (PSF) of ACS is very stable in time. However, 
coma, astigmatism and defocus from surface height vari- 
ations of the two CCDs lead to variations over the field of 
view that need to be taken into account. The variations 
are much weaker than for the WFPC2 camera but remain 
non-negligible. Also the 'breathing' of HST changes the 
focal length which leads also to a small time variability in 
the PSF. In Figure El we show the mean ACS PSF com- 
piled from ^^500 stars in the extensive gems area and 
the variation of the PSF over the FOV. 

The variations were also investigated bv lKrist fc Hookl 
(|2003i) from the crowded field of 47 Tucanae. Such an 
analysis was not possible for individual gems fields due 
to the small number of stars (~10) per field. Thus a si- 
multaneous characterisation of the spatial and temporal 
variations was not possible. However, since the spatial 
variations dominate (further details of our investigations 
of the ACS PSF variations will be given in a dedicated 
technical paper, Jahnke et al. 2004, in prep.) we used the 
large number of unsaturated stars in the gems area to 
construct an empirical PSF individually for each object. 
At a given position we combined the nearest ^^35 undis- 
turbed stars to create a position-specific PSF estimate. 
In this way we average over time, but only stars from 
a radius < 40" were used, and PSF shape errors due to 
spatial variation were minimized, while a very high S/N 
was achieved for each of these PSF estimates. Finally, 
the subimage of each AGN and its connected PSF were 
registered to a common centroid. 

From our PSF analysis we found that while coherent 
large-scale variations were essentially absent within each 
stack of '--^35 PSF stars, there was still considerable mis- 
match between the individual stars, in particular in the 
central pixel regions. As such mismatched pixels could be 
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TABLE 1 
Objects of the sample. 
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''ID from the COMBO- 17 catalogue 
''GEMS tile number (1-63 GEMS, 80-95 GOODS region) 
"^Photometric redshift from COMBO-17; for COMBO-17 19965 a re- 
vised redshift from an ongoing objective prism survey was used 
''ij-band magnitude (Vega zeropoint) from COMBO-17 
•^Total F606W- and F850LP-band magnitudes (AB zeropoint) as 
measured with the ACS. Photometric errors range from 0.07-0.10 
for both filter. Three objects are imaged on two different tiles. 
Due to nuclear variability this can result in different F606W and 
F850LP-band fluxes. 



spuriously assigned to a host galaxy, we took the vari- 
ations within each PSF stack to derive rms frames de- 
scribing an inherent PSF uncertainty; these were then 
also included in the variance images, artificially reducing 
the weight in the inner pixel regions. PSF rms errors per 
pixel range from up to 30% in individual pixels inside 
1 FWHM to 5-15 % inside 0'.'2, and to generally below 
5% outside. 

3.3. Peak scaled PSF subtraction 

For luminous AGN at high redshifts, separating the 
galaxy image from the nuclear point source is a daunting 
task. Even with a very good knowledge of the PSF, the 
problem is still that the relative scalings of galaxy and 
AGN are not known a priori. In fact, one cannot even 
be certain that the host galaxy is detectable at all and 
not swamped by the central point source. We therefore 
started out with the well-established technique of simple 
PSF subtraction. 

In each case, the PSF was scaled to the central flux 
of the AGN, integrated inside a circular aperture of 4 
pixel (0'.'12) diameter centered on the nucleus. This ra- 
dius encircles approximately 34% of the total energy of 
a point source. A smaller radius (e.g. 1 pixel) would 
become too sensitive to shot noise and PSF mismatch. 



while larger radii contain a significantly higher fraction 
of the total flux (50%/60%/70% at 3/4/5 pixel radius) 
and thus would make a detection of any host galaxy com- 
ponent successively harder. 

This procedure somewhat oversubtracts the nuclear 
component by an amount corresponding to the under- 
lying host galaxy contribution inside the encircled re- 
gion. However, this method yields strictly conservative 
estimates the host galaxy flux, i.e. always unserestimat- 
ing it. We used extensive simulations to determine cor- 
rection factors for this oversubtraction (see Section IX^ . 
The peak-scaled PSF subtraction method has the advan- 
tage of being independent of any assumption about the 
host galaxy morphology. 

3.4. Two-dimensional deblending 

As a second method we e mployed the mod elling pack- 
age GALFIT in Version 1.7a l)Peng et alJl2002|) that allows 
the simultaneous fitting of several two-dimensional com- 
ponents to an image, convolved with a given PSF. We 
describe our application of GALFIT to quasar images in 
Sanchez ct al. (2004) in detail. 

In total we ran GALFIT in three configurations, always 
fitting two components, which were t he point-source nu- 
cleus and either an exponential disk llFreemanlll970l) . a 
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Fig. 2. — Top: PSF in the F606W-barid, shown with three display levels to enhance structures (size is 3'.' 84x3'.' 84). Bottom: 
Subtraction of a total mean PSF from the PSFs in each of 4x4 subtile regions to show the spatial variations of the PSF in the 
inner region (size is l'.'92x l'.'92). 



de Vauc ouleurs spheroid ijde VaucouleurEll948() . or a 
ISersid l|1968f) model with the Sersic index as a free param- 
eter. However, simulations (Sect. I3.6|l showed that with 
the present data, the GALFIT version used^ was operat- 
ing near its limit, due to the very high contrast between 
nuclei and host galaxies. In the context of this paper 
we thus used GALFIT only as a cross-check on the peak 
subtraction method. 

3.5. Detection sensitivity 

To determine the limits for detecting host galaxies we 
constructed a sample of 200 randomly selected unsatu- 
rated stars, 100 in each observed band, to mimic unre- 
solved quasars. This way we could investigate how our 
nucleus-removal techniques responded to a undetectable 
host galaxy, and we could set limits on the size and shape 

® The GALFIT version 1.7a was recentering the PSF to given 
coordinates using a convolution with a narrow gaussian, not by 
shifting the PSF by means of rebinning to a new position. The 
latter is better for the application to AGN decomposition and is 
now incorporated in later versions of GALFIT. 



of expected residuals, thus lower flux limits for detectable 
host galaxies. For these stars the PSFs were created in 
exactly the same way as for the AGNs. The object itself 
was always excluded in the PSF production thus each 
test star and its PSF were fully independent. 

This set of simulated 'naked quasars' showed that in 
88% (97%) of all cases, any residuals - which could be 
taken as spurious (g)host galaxy detections - had fluxes 
of less than 5% (10%) of the total object flux. From this 
we adopted the condition that a real detection should 
show a residual flux after peak subtraction of at least 5% 
of the total flux, corresponding to a maximum nucleus- 
to-host ratio of 20. Because of the systematic oversub- 
traction inherent in the procedure corrections for the flux 
have to be applied (see next section). 

The final decision if a host galaxy is resolved is based on 
this criterion. In addition we visually inspected whether 
the detected flux indeed came from a host galaxy or 
whether other, unmasked structures were present, using 
the peak subtracted images and radial profile. If this 
could be ruled out we classified a host galaxy as detected. 
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Fig. 3. — Simulation results on systematic offsets and errors in recovered host galaxy magnitude. Left: Recovered host magnitude as 
function of input magnitude for a range of morphologies, differing in scales and nucleus-to-host flux ratios matching our sample, for the 
two bands. The adopted reliability regions end at the dotted straight lines (see text). The F606W- and F850LP-bands are simply shifted 
with respect to each other by their zero-point offset (ZPpggew ~ ZPpgsoLP = 1.59). Right: Difference of in- and output magnitudes Am 
vs. output F606W band magnitude of the mean relation to the left. The overplotted errors are the la error ranges of the measured spread 
of host galaxy flux, as determined from the set of simulations with different parameter combinations and different noise realisations (see 
text). The vertical lines gives the edge of the reliability region for the F606W band as above. For the F850LP band again the zeropoint 
has to be shifted. 



3.6. Systematic offsets and errors 

While the mere detection of an AGN host galaxy can 
be achieved with comparably little effort, the determina- 
tion of flux error bars and systematic offsets is much more 
complicated. We performed extensive simulations of arti- 
ficial quasar images, composed from empirical PSFs and 
host galaxy models plus artificial noise matching the ac- 
tual flux and noise distribution in real images. These sim- 
ulations are described in detail bv lSanchez et alJ |2004). 
We applied the peak subtraction nucleus removal as well 
as GALFIT to a set of ^^2000 quasar images created in 
this way. Comparing input and output parameter val- 
ues yielded mean magnitude offsets as well as statistical 
errors for the individual host galaxy magnitudes (Fig. ^ . 

The simulations give reliability regions and error bars. 
The left panel in Figure 13 shows which magnitudes are 
recovered for a given synthetic host galaxy. Since the 
input set covers a large range of different morphological 
configurations, scale lengths, nucleus-to-host ratios, etc., 
the recovered values will scatter. Close to the detection 
limit, the scatter and the corrections grow rapidly as a 
function of magnitude; additionally, the ability to differ- 
entiate between different morphological types will gener- 
ally be lost. The combination of these effects is refiected 
in the spread of the output of the simulations. This mea- 
sured spread is a direct estimate for the uncertainties of 
the total flux (right panel in Figure 

From these simulations we adopt approximate regions 
in brightness where host galaxy magnitudes can be reli- 
ably determined, with correction of 0.25 to a maximum 
of 0.6 mag. These regions go down to F606W = 26.2, 
F850LP = 24.6 for the peak subtraction method and 
the present data. Outside the corrections and errors in- 
crease. In three cases, marked in Table |5] with '?' in 



column Ziig, the observed magnitudes extend to outside 
these regions; here we continue using the derived cor- 
rections for these three objects, but the so derived host 
galaxies are more uncertain and their magnitudes should 
be taken with care. Notice that the F850LP band data 
are substantially shallower than the F606W band, mainly 
a consequence of the ACS detector sensitivity. 

4. RESULTS 
4.1. Detected host galaxies 

Using the above criteria for detecting a residual host 
galaxy, we find nine of the 23 host galaxies to be resolved 
in both bands, although some lie close to the sensitivity 
limit. One object formally fell above the 5% level in one 
band but not the other; COiviBO-17 33630 at z = 2.719 
might be marginally resolved in F850LP and shows a 
structure at l"distance that might be a tidal arm or a 
foreground object. With this object lying at the highest 
redshift of the sample we do not consider this a clear 
detection. As mentioned, for three further cases in the 
F850LP band the host galaxies are very faint (marked 
with a '?' in the Zhg column in Tab.[2Jl. While their fiux 
is above 5% of the total, their raw magnitudes fell 0.8- 
1.0 mag outside the reliability region where corrections 
and associated errors are still small. This low S/N is also 
reflected in the radial profiles (see Appendix 0. 

As described above, tests with field stars show that 
12% of all objects (^--^3 objects) might show spurious 'host 
galaxies' at the 5% flux level, and 3% (0 or 1 objects) at 
the 10% flux level. In F606W five of our objects faU 
with their host fiuxes between these two values. In the 
F850LP-band these arc four which include the three un- 
certain ones from above. According to statistics 1-3 of 
these might be spurious detections. However, including 
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Fig. 4. — Peak subtraction host galaxy magnitudes. Top: uncorrected for oversubtraction, bottom; corrected. Left F606W, 
right F850LP. The arrows give upper hmits for the unresolved HGs, the open symbol for F606W is the HG extracted from the 
stacked unresolved AGN images. The horizontal error bar for the stack gives the la redshift range of the coadded AGN. In the 
bottom diagrams the vertical errors are uncertainties from the subtraction of the nucleus and are given as determined from the 
simulations (see Sect. IH.fil and Fig.|3l. 



or excluding these more uncertain data points in the fol- 
lowing analysis does not have an influence on the conclu- 
sions drawn. 

For each object the host galaxy flux was determined 
by simple aperture photometry after subtraction of the 
scaled PSF, excluding resolved companion objects. The 
radius of the aperture was matched the used image size. 

All extracted magnitudes are collected in Table El and 
shown together with the upper limit for the unresolved 
objects in the top panels of Figure01 The extracted host 
galaxy images and radial surface brightness profiles are 
shown in Fig. lAllI in the appendix. To illustrate the 
behaviour of true point sources, we included a selection 



of 24 hel d star s, 12 each in the F606W- and F850LP-band 
fFigures |A12l and lA13p . that were subjected to the same 
PSF determination and peak subtraction as described in 
Section IX5l We plot the same profiles as for the AGN. 
This selection is random apart from the fact that the 24 
stars were observed on 24 different tiles. For most of the 
stars there is no systematic positive residual flux visible, 
as expected. The few that do show positive fluxes form 
the spurious detection statistic described above. 

While with deep er images or at lower redshift 
ijSanchez et alJl200^ the morphological appearance of 
the host galaxies can be determined, this was generally 
not possible for the present data. Apart from the ex- 



8 



K. Jahnke et al. 



Composite 'stacked' object, z ~ 2.3 
original F606W nuc. subtracted 




r (arcsec) 



Fig. 5. — The stacked image of individually unresolved host 
galaxies. Top: grey scale original image and PSF subtracted host 
galaxy image (both show 2" X 2") in the F606W-band. On top of 
the linear grey scale plot logarithmic contours with 0.5 dex spacing 
are overplotted. Bottom: radial profiles of F606W-band. The 
curves show the upper data points with error bars from the original 
image, the PSF (dashed red line) and peak subtracted host profile 
(lower points with error bars from bootstrapping). The host galaxy 
is shown without correction for any oversubtraction of the nucleus. 



tentions in COMBO- 17 00784, apparent residual struc- 
ture visible in the colour images shown in the Appendix 
is dominated by the PSF subtraction procedure. This 
includes elongations and apparent ofF-centering that is 
due to limited centering precision of the order of 0.05- 

0. 1 pixel (see e.g. combo- 17 05696 from tile 2). Thus 
we adopted the host galaxy morphological class in all 
cases as 'undecided', and then apphed the systematic 
corrections for oversubtraction provided by the simula- 
tions (Fig.OJ. The corrections were typically of the order 
of ~ 0.4 mag, as documented by the columns with 'cor' 
subscripts in Tabled Here we also list the estimated un- 
certainties resulting from our extensive simulations. The 
distribution of corrected magnitudes is displayed in the 
lower panels of Figure^ 

4.2. Detecting the undetected hosts 

For 14 AGNs we find that the individual PSF- 
subtracted residuals are consistent with non-detections, 

1. e. the magnitudes of individual host galaxies lie below 
the 5% limit. These objects are marked by arrow sym- 
bols in Figures 0] and |H1 It is interesting to note that 
these flux limits on the hosts are by no means all out- 
standingly faint. For many, the reason why they were 
not detected is the high contrast with the AGN which 
for these objects has above average brightness. Only 
three objects show non-detections that indicate excep- 
tionally faint host galaxies. We will further discuss the 



implications of the detection limits in Section l5.1l 

In an attempt to assess at least the mean host galaxy 
properties of the unresolved AGN, we simply coadded 
the images of all 14 objects, one of which was observed 
in two frames. This yielded a very deep image with 
effectively 15 orbits of integration time. We also com- 
bined their PSFs, weighted by the relative flux of the 
AGNs, and created combined variance and PSF rms 
frames. The higher S/N in combination with the lower 
PSF noise (PSF and AGN position are sampled at 14 
different subpixel points) yields a nominal increase in 
surface brightness sensitivity of ^^1.5 mag. The resulting 
'object' has a mean redshift of z = 2.3 (weighted by AGN 
flux), and total magnitudes of mpeoew.tot = 21.5 and 
?7iF850LP tot = 21.6, respectively. This image is shown in 
Figureinj With the higher sensitivity we now indeed find 
a host galaxy component in the F606W-band image after 
PSF subtraction of 4.4% of the total flux. The radial sur- 
face brightness profile also shows a small excess over an 
unresolved point source. In both cases this fiux is highly 
significant as we confirmed using a bootstrap simulation 
for the composition of the coadded frame from the 15 
frames. In the bootstrap simulation we constructed 100 
new sets of 15 frames each, drawn with repetition from 
the original 15 frames, coadded the images in each set 
and did the flux analysis as above. The uncertainty in 
the total flux estimated from these 100 realisations is 
(7 = 1.05% of the total flux, or 25% in host galaxy flux. 
All realisations yielded substantial positive fluxes. The 
error is resulting from a combination of PSF uncertainty 
and the noise inside the scaling aperture of 4 pixel di- 
ameter. We show the uncertainties in the radial sur- 
face brightness determined from bootstrapping as error 
bars for the derived host galaxy in Figure El The so ex- 
tracted magnitudes for host galaxy and nucleus in the 
F606W-band are listed in the last row of Table 12 (the 
'stack' object). The F850LP-band stack, however, with 
its lower sensitivity showed a much weaker signal than 
the F606W-band, too faint to reliably be classifled as 
resolved. 

5. DISCUSSION 
5.1. UV colours 

From the corrected magnitudes we have derived 
(F606W - F850LP) colours for the detected host galax- 
ies. These are listed in Table |21 now including also 
correction for Galactic dus t extinction. However, with 
E{B -V) ^ 0.008 ( Schlege ]"eral]ll998|) these values of 
^(F606W) = 0.024 and A(F850LP) = 0.014 negligibly 
affect the colours. 

At these redshifts, the observed photometric bands 
correspond to the rest frame ultraviolet, ranging from 
2160 A at z = 1.8 to 1616 A at z = 2.75 in F606W and 
3150 A to 2350 A in F850LP, respectively, so for this red- 
shift range pure rest-frame UV colours are observed. Fig- 
urel^lshows the measured values plotted against z. There 
is no discernible colour trend with redshift. All points 
fall within a relatively narrow range of colours; apart 
from COMBO-17 51835, the objects occupy a band of 
-0.2 < (F606W-F850LP)obscrvcd < 1-0. The open sym- 
bol represents our stacked 'average' AGN constructed 
from the 14 unresolved objects. Although it was not re- 
solved in the F850LP band, the upper limit on its colour 
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TABLE 2 

Photometry results for the F606W- (for brevity: 



V"-) AND F850LP- {'Z') BANDS. 



ID 


Tile 


Vtot" 


ng 


ng,cor 






Zt t° 




hg.cor 




N/Hz.cor" 


19965 


23 


20.59 


23.9 


23.3±0.05 


20.7 


11.3 


20.29 


23.8 


23.5 ±0.2 


20.4 


17.5 


30792 


82 


21.60 


24.6 


24.4±0.2 


21.7 


12.0 


22.06 


24.2 


23.9 ±0.2 


22.3 


4.2 


18324 


19 


22.00 


23.6 


23.0±0.05 


22.5 


1.6 


21.33 


22.7 


22.4 ±0.2 


21.9 


1.6 


05498 


01 


22.91 


25.5 


25.2±0.2 


23.1 


7.1 


22.28 


25.0? 


24.4 ±1.2 


22.5 


6.0 


51835 


55 


23.01 


26.0 


25.7±0.3 


23.1 


11.2 


22.41 


23.8 


23.4 ±0.2 


23.0 


1.6 


00784 


05 


23.28 


25.0 


24.7±0.2 


23.6 


2.8 


22.80 


24.0 


23.7 ±0.2 


23.4 


1.2 


05696 


02 


22.73 


24.7 


24.4±0.2 


23.0 


3.6 


22.32 


23.8 


23.5 ±0.2 


22.8 


1.9 


05696 


03 


23.14 


24.2 


23.9±0.15 


23.9 


0.9 


22.68 


24.0 


23.7 ±0.2 


23.2 


1.6 


07671 


07 


22.35 


25.6 


25.2±0.2 


22.4 


13.2 


22.24 


25.0? 


24.4 ±1.2 


22.4 


6.3 


07671 


15 


22.36 


25.5 


25.1±0.2 


22.5 


11.9 


22.22 


24.6 


24.3 ±0.3 


22.4 


5.5 


11922 


11 


22.65 


24.9 


24.6±0.2 


22.9 


4.8 


21.93 


25.2? 


24.5 ±1.5 


22.0 


9.3 


stack 




21.47 


24.9 


24.3±0.05 


21.6 


12.3 


21.60 


unresolved 







Note. — All magnitudes are uncorrected for Galactic extinction, A(F606W) = 0.024 and A(F850LP) = 0.014. 



^Total magnitude of the object. 

''Raw measured host galaxy magnitude. 

'^Host galaxy magnitude corrected for oversubtraction. 

"^Corrected nuclear magnitude. 

•^Nuclear to host galaxy flux ratio. 

'objects marked '?' lie outside the adopted reliability region (see 
text) . 



TABLE 3 

Observed (F606W — F850LP) colours and star formation 

RATES FOR THE RESOLVED HOST GALAXIES AND THE 'STACKED 
AGN', IN CASE OF THE COLOURS BEING PRODUCED BY CONSTANT 
STAR FORMATION. 



ID 


Tile 


2 


(F606W - F850LP)'' SFR F606W 
[Mo/year]b 


19965 


23 


1.90 


-0.2±0.2 


11 


30792 


82 


1.929 


0.5±0.3 


4 


18324 


19 


1.990 


0.6±0.2 


15 


05498 


01 


2.075 


0.8±1.2 


2 


51835 


55 


2.179 


2.3±0.4 


1.5 


00784 


05 


2.282 


1.0±0.3 


4 


05696 


02 


2.386 


0.9±0.3 


6 


05696 


03 


2.386 


0.1±0.3 


9 


07671 


07 


2.436 


0.8±1.2 


3 


07671 


15 


2.436 


0.9±0.4 


3 


11922 


11 


2.539 


0.1±1.5 


5 


stack 




2.3 


< 0.0 
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100% z=5 formation 



0.1% lOOMyr 



1% 



10% 



i 



100% 
lOOMyr 



''Galactic extinction corrected. 

'^Star formation rates inferred from the UV flux in the observed 
F606W-band. 



is actually consistent with the values derived for several 
of the detected objects. 

One critical issue in measuring UV luminosities of 
barely resolved AGN hosts is the lingering possibility 
that flux may have spilled over from the nuclei. This 
could have happened as a purely observational artefact 
due to imperfect PSF removal, or physically by scattering 
of UV photons off dust in the the host galaxy. The lat- 
ter phenomeno n is known to be rele vant in high redshift 
radio galaxies ijVernet et al.l I200H) . Independently of 
the underlying mechanism, any such cross-contamination 
should be visible in a correlation of host galaxy with 
nuclear colours. Figure [7\ shows these colours plotted 
against each other. No correlation is visible and we con- 
clude that a substantial contamination of the host galaxy 



1.8 2 2.2 2.4 2.6 

redshift z 

Fig. 6.— Observed colours (F606W - F850LP) of the sam- 
ple from PSF peak subtraction (circles), the open symbol marks 
the upper limit for the 'stacked' AGN. Overplotted are tw o sin- 
gle burst models from (solar metallicitv IBruzual Hz Chariot 2003) 
(solid lines). The upper curve is for a passively evolving burst at 
z = 5, the lower for burst of 100 Myr age, relative to each redshift. 
The dot-dashed lines are mixtures between the two, with a (from 
top) 0.1%, 1% and 10% fraction of mass of the 100 Myr population 
on top of 99.9%, 99% and 90% of the 2 = 5 population. 



light from the AGN is very unlikely. Notice also that 
when considering physical scattering in the hosts, pow- 
erful radio galaxies are huge massive entities of generally 
vastly different appearance compared to the relatively 
modest AGN hosts featuring in our sample. 

UV colours can be converted to UV spectral slope /? 
independent of redshift, when assuming that the SED 
can be described in the form F\ oc A*^. With (3 known. 
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Fig. 7. — Nuclear vs. host galaxy colours for the objects in 
the sample with available colours. The two objects imaged twice 
are connected with a dotted line (in one case not visible to near 
identical colours). The dashed line is the 1:1 relation to guide the 
eye. 



the absolute magnitude at 200 nm, M2oonm, can be com- 
puted directly from the F606W- or F850LP-band appar- 
ent magnitudes. We do this for both the host and the 
nucleus and these values are shown for the full sample in- 
cluding upper limits in Figure|Sl which illustrates the two 
main constraints for resolving a host galaxy, apart from 
compactness (see Sect. l4.2|) . The reliability of host galaxy 
photometry is constrained by S/N, thus dependent on 
the filter which have different depths. This is marked 
by the horizontal dashed lines which show the M2oonm 
magnitude of host galaxies at the F606W = 26.2 mag 
and F850LP — 24.6 mag edges of the adopted regions of 
reliability (see Section l3.6|l . at mean redshift and mean 
/3. 

As the second effect the maximum nucleus-host con- 
trast appears as the scatter of the unresolved objects 
around a line shifted by 3.2 mag from unity, correspond- 
ing to 5% of the total flux. Here the scatter is only in- 
duced by the assumption that all unresolved host galax- 
ies have a the /3 value the upper limit determined for the 
stacked image. The diagonal dotted lines mark lines of 
10%, 20%, 33% and 50% of the total flux associated with 
the host galaxy. Thus in total the region right of the solid 
diagonal line is inaccessible to host galaxy detection with 
the current data and method of analysis. We would like 
to emphasise that the similarity of nuclear properties for 
resolved and unresolved host galaxies suggests that also 
the host galaxy properties are similar - thus supporting 
that the data point for the coadded stacked object is not 
far off the individually resolved objects in all plots. 

In Figure El spectral slope is plotted against M2oonm 
for the host galaxies, compared to the mean value for a 
sample of 794 Lyma n break galaxies at redshift z ^ Z 
l|ShaDlev et al.ll2003D . all values uncorrected for the in- 
fluence of dust in the galaxies. 

The slight anticorrelation that seems to be visible be- 
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Fig. 8. — Nuclear vs. host galaxy absolute magnitude at 200 nm. 
Shows are objects with resolved host galaxies (solid symbols), the 
coadded stacked object (open circle with upper limit arrow) and 
upper limits for objects with unresolved host (arrows). The diago- 
nal lines mark positions of constant fraction of host galaxy light of 
5%, 10%, 20%, 33% and 50% of the total light, assuming a constant 
spectral slope /3. The horizontal dashed lines show the magnitude 
of host galaxies at the F606W = 26.2 mag and F850LP = 24.6 mag 
detection limit, at mean redshift and mean /3. 



tween /3 and UV luminosity suggest that the more lumi- 
nous host galaxies are having the steeper spectral slopes 
or bluer colours. This would mean that more luminous 
host galaxies had bluer colours and thus more UV light 
from young stars. To test this we computed the Spear- 
man rank-order coefficient for this data set. The test 
gave a probability for the zero hypothesis - uncorrelated 
data - of 9%. Thus the zero hypothesis can not even be 
rejected on a 2cr level and thus the anticorrelation is not 
significant. 

A few of the most UV-luminous host galaxies fall into 
the same magnitude-colour regime as the LBGs, while a 
number of objects are substantially redder. Thus we do 
not find a positive correlation of the amount of host UV 
light and luminosity. 

5.2. Origin of the host galaxy UV flux 

The top solid line in Figure El shows the theoretical 
colour of a galaxy that formed all its stars at z = 5 
and evolved passively afterwards. This colour was com- 
pute d from single stellar popu lation (SSP) models taken 
from lBruzual &: Charloti l|2003>l . The assignment of so- 
lar metallicity in this context is arbitrary and subject 
to discussion (see below). However, qualitatively there 
is no strong dependence on metallicity, the exact forma- 
tion redshift o r the particular c hoice of IM F for a given 
model family IjSalneteil IT955I or iChabrieri 1200 3) in this 
wavelength range. 

This comparison shows that the measured UV colours 
of our detected AGN hosts as well as of the stacked 
'mean' host galaxy are markedly bluer than expected for 
an 'old' population at that epoch, i.e. tage ^Hubble- Any 
correction for dust in the host galaxies will strengthen 
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Fig. 9.— Spectral slope 13 in the UV {Fx oc A'^) vs. 200 nm 
absolute magnitude, uncorrected for dust. The open symbol marks 
th e 'stacked' object, the crossed square is the mean value obtained 
bv lShaplev et aEl 120031) for 794 Lyman break galaxies at z ~ 3. 



this result. Thus the blue light detected in these galax- 
ies must come from relatively young stars. These stars 
could be forming continuously, in which case the UV lu- 
minosities can be interpreted as indicators of the star 
formation rate in these galaxies. We follow up on this 
option in the next subsection. Alternatively, the UV flux 
could be the afterglow of a past starburst. This is the 
option we consider first. 

5.2.1. Recent starburst 

Given the unknown dust absorption and metallicities 
in the host galaxies, a single UV colour is insufficient for 
performing a detailed age dating of starburst or a de- 
composition of stellar populations. However, we want to 
compare to the available theoretical colour range spanned 
by galaxy formation at z = 5 and a very recent (100 Myr) 
starburst to illustrate how mixing of an old underlying 
population with most of the mass and a recent starburst 
influences the colour. In Figure the solid lines mark 
these extremes. In between these the dot-dashed lines 
show how contribution of a 100 Myr component would 
influence the colour of an otherwise old population. From 
top to bottom 0.1%, 1% and 10% in mass are added to 
the old population. There is a strong degeneracy be- 
tween the choice of burst age and the mixing ratio. For 
10 Myr less than one tenth in mass is required to pro- 
duce the same UV colour compared to 100 Myr. If we 
choose 100 Myr as a timescale similar to the dynamical 
timescale in galaxies, the masses involved in that star- 
burst would be of the order of a few percent of the total 
stellar mass. 

For the assumption of only one single-aged population 
we can rule out a very high formation redshift - z = 5 
corresponds to ages of 3.5-2.5 Gyr a,t z — 1.8-2.5 in 
the chosen cosmology. For the adopted set of models, 
the resulting age estimates would range mostly between 



-0.1 and -0.7 Gyr. 

We note that the colour tracks in Figure El are largely 
flat over the redshift range of interest, and that the two 
pure population and the different mixing ratios corre- 
spond to different colours almost independently of z. 
Since our measured colours are all quite similar, we con- 
clude that the luminosity weighted ages of the the UV- 
dominating stellar population must be rather similar, un- 
less younger ages and more reddening in some objects 
conspire. 

Clearly, a single UV colour is insufhcient to perform 
a reliable age dating, with all broad band colours be- 
ing affected by various degeneracies with respect to dust 
and metallicities. However, we have reason to believe 
that at least the central lines of sight towards the AGN 
are reasonably free of dust extinction (because the AGN 
sample is selected by optical/UV flux), and we therefore 
do not expect dust to play a major role. At any rate, 
significant dust extinction would make the host galaxies 
intrinsically bluer than what we observe. On the other 
hand, assuming a metallicity lower than solar would shift 
all curves in Figure |H| downward, resulting in older age 
estimates. Reducing the metallicity to Z = 0.004 (1/5 
solar) gives an single burst age increase by a factor of 
two (for Bruzual & Chariot models). We conclude that 
if the UV light in our host galaxies is emitted by a pas- 
sively evolving population of young stars, this population 
is typically much younger than a Gyr. 

The diagram in Figure |H1 shows that no correlation ex- 
ists between the amount of stellar UV light from the host 
galaxies and the amounts of UV radiation produced by 
the nuclei. If the latter is taken as a measure of the 
amount of matter accreted by the nuclei then in the con- 
text of a recent starburst the size of the starburst and the 
amount of accreted matter must be governed by different 
mechanisms. If the accretion rate is primarily defined by 
the nuclear mass, and a correlation between galaxy and 
black hole mass is assumed, then the size of the star- 
burst is independent of the bulk stellar mass of the host 
galaxy. Other factors must be dominating the amount 
of (gas) mass involved in the starburst. This can be ei- 
ther the total amount of gas available in the galaxy, the 
size of the region involved in the starburst, the strength 
of an interaction responsible for the starburst, etc. In 
any case the amounts of gas involved are variable for 
a given nuclear luminosity. Including the upper limits 
in Figure |H1 the host luminosities span e.g. —4 mag at 
-^^200 nm(nuc) ~ —22, i.e. the amount of gas involved 
can vary by a factor of — 40 or more. 

5.2.2. Estimating a host star formation rate 

We now interpret the detected UV emission in the al- 
ternative framework of being due to young stars forming 
continuously in the AGN host galaxies. Under this as- 
sumption it is possible to estimate the star formation 
rate (SFR) of the host galaxies fr om the measured rest 
frame UV luminosities. Following iKennicuttI ijlQQST ) and 
using the conversion of AB magnitudes into monochro- 
matic fluxes we obtain 



SFR 



Mo 
year 



= 1.8 X 10 



-27 



^2-i^Q-0.4(mAB+48.6) 



where di is the luminosity distance to the AGN in cm, 
and mAB is the observed UV magnitude at an arbitrary 
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Lyman break galaxies 



i 



1.8 2 3.2 2.4 2.6 

redshift z 

Fig. 10. — Rest frame 200nm luminosities, and star formation 
rates as derived from from the F606W-band, both uncorrected for 
dust. The open symbol marks the SFR of the 'stacked' object 
created from the AGNs with individually unresolved h ost galax- 
ies. T he horizontal dashed line is the value obtained bv lErb et al.l 
120031) for Lyman break galaxies at z = 2.5. 



wavelength between 1500-2800 A. With this formula we 
can now convert our F606W band luminosities (which 
are much deeper than the F850LP band data) into star 
formation rates. As long as intrinsic dust attenuation is 
neglected, these values are of course mere lower limits. 

The resulting SFR values are listed in Table 13 and 
generally amount to a few solar masses per year, with 
remarkably little variation. The maximum value found 
is ~15 Mq yr~^, and the mean is ~6 Mq yr~^, including 
the results from the stacked image of individually unre- 
solved host galaxies. This value is thus representative for 
the full sample of 23 AGN. We show the distribution of 
star formation rates vs. redshifts in Figure [TOl No trend 
emerges, consistent with the previously established ob- 
servations that neither colours nor luminosities show any 
significant trend with z. 

This is again compared to the star formation rates for 
LBGs at z = 2 ... 2.6 from Erb et al. ( 2003), as deter- 
mined from the UV flux and uncorrected for dust. There 
are a few host galaxies with higher UV flux while for the 
majority it is smaller by a factor of 2-3. 

However, there is a principal caveat in our guiding as- 
sumption of this subsection: The UV light from galaxies 
with strong ongoing star formation is expected to be to- 
tally dominated by the youngest stars. lKennicufr3 l|1998fl 
pointed out that the spectral shape of galaxies with a 
constant SFR over at least ~ 100 Myr is basically flat 
(in /,,) between 1500-2800 A, assuming a SalDcter ( 1953[) 
IMF. This would lead to an expected UV colour for our 
objects of (mpeoew — ?tif850Lp) = or spectral index 
P — —2, more or less independently of redshift (within 
the z range of our sample) , inconsistent with our observa- 
tions for most objects (see previous subsection. Figure El 
and Table OJ. In other words, most detected host galax- 
ies of our sample have UV colours that are too 'red' for 
a simple continuous star formation scenario, while the 



stacked host galaxy is roughly consistent. If for the re- 
solved hosts the F850LP band fluxes were used instead of 
the F606W band, SFR values would be higher by roughly 
a factor of 2. 

This apparent inconsistency could be resolved in sev- 
eral ways. The initial assumption could be wrong, and 
the UV flux originates not from freshly formed stars, 
but from a passively evolved starburst as outlined in 
Sect. 15.2m The star formation rate might not have been 
constant over the past, so that varying amounts of stars 
of different masses and ages would have been formed; 
such a configuration is always possible, and our only ar- 
gument against it would invoke Occam's razor. Finally, 
as discussed in the previous section, there could be an 
underlying older stellar population contributing more to 
F850LP and less to F606W (but note that by sample de- 
sign, also F850LP is completely below the Balmer jump 
at all relevant redshifts). This is the scenario that is 
favoured by our inte rmediate redshift data from gems 
(iSan chez et !liI] l2004D . 

Clearly, in this study our current data set of just two 
UV bands is insufficient to settle this ambiguity. How- 
ever, in all three cases there would be no continuous star 
formation as inferred for LBGs. The blue light would be 
a result not of a continuous process but of one or more 
events in the past of the host galaxy that triggered star 
formation. Whether galaxy interaction or merger inci- 
dents were responsible or the formation of bars or spiral 
arms is involved can not be investigated with the present 
set of data. 

5.3. Comparison with other AGN host galaxy studies 

Even at low redshifts, colour data of AGN host galax- 
ies are relatively scant, except for rather low- luminosity 
AGN where the host galaxy can be separated with 
relative ease. In those cases, colours were gener- 
ally found to be consistent with morphological types, 
in particular for the prevailing disk-type ho st galaxies 
l)Kotilainen fc WarJl994HSchade et"ani2000() . However, 
when higher nuclear luminosities are observed (which 
generally correlate with a larger bulge component), 
then a tendency towards abnormally blue colours and 
younger stellar populat ions emerges ijKauffmann et af] 
l2003H Jahnke et alJ200J) . This tendency appears to hold 
also at intermediate redshifts, as demonstrated in our 
companion HST paper (Sanchez et al. 2004), where we 
find more than half of the investigated host galaxies to 
have bluer colours than what would be expected from 
their morphological types. However, we can not confirm 
increasing amounts of blue light from young stars with 
incre asing host luminosity as found by ( Kauffmann et aQ 
I2OOI . Our data are consistent with constant UV flux for 
all hosts. 

At high redshifts (z > 2), colour information is avail- 
able only for a handful of objects. In a ground-based 
study of six bright radio-loud quasars, Lchnert et al} 
(1992) found indications that the hosts were very blue, 
actively star-forming galaxies. For three high lumi- 
nosity quasars , one radio- loud and two radio-quiet, 
lAretxaga et alJ (1199 8*) claimed very luminous envelopes 
and star formation rates of se veral hundreds s olar m asses 
per year. More recently, Hutchin gs et al.l (|20O2') pre- 
sented optical HST observations of three radio-loud and 
two radio-quiet quasars and found less extreme, but qual- 
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itatively similar results. 

Because of their high radio power and optical luminosi- 
ties, many of the quasars observed in previous studies are 
probably quite different from the moderate-luminosity 
hosts of moderate-luminosity radio-quiet AGN that we 
focus on. It is nevertheless interesting to see that the 
presence of a significantly enhanced UV continuum and 
young stars seems to be at least qualitatively similar be- 
tween QSOs of intermediate and high luminosities. 

6. CONCLUSIONS 

We performed the hitherto largest study of host galax- 
ies properties of a complete sample of high-redshift AGN. 
We detected the hosts and extract colour information in 
9 of the 23 AGN, and we also achieved a statistical detec- 
tion of the host in the remaining 14 from a stack analysis. 
The UV luminosities can be interpreted in three ways: ei- 
ther as contribution from a passively evolving population 
of relatively young stars, forming typically 0.5 Gyrs ago, 
as a mix between a population of old (e.g. tage ~ inubbic) 
stars and a small contribution of a recently formed young 
population (e.g. 0.1%-10% in mass at an age of 100 Myrs 
or 1/lOth of this for age 10 Myrs), or as an indicator of 
ongoing star formation at a level of ~2-15M0 yr~^ (un- 
corrected for internal dust attenuation). While the first 
possibility is very simplistic and appears unphysical, the 
UV colours actually favour the two burst interpretations; 
but the possibility of on-going star formation cannot be 
completely ruled out from our data. 

In the framework of combined old and young pop- 
ulations, it is remarkable how similar the host galaxy 
colours are within the sample, and, unless different mass- 
age combinations conspire, hence the estimated stellar 
mass fractions and ages. The derived young popula- 
tion mass fractions and ages are also very similar to 
the values estin iated in our compa nion gems study of 
AGN at z < 1 l|Sanchez et al.ll2004j) . where we find ab- 
normally blue rest-frame U — V colours for a substan- 
tial fraction of host galaxies, particularly the most lu- 
minous AGN in the sample. Even more, these colours 
and ages are in turn very close to the mean values ob- 
tained from our gro und-based low-z multicolour sample 
l),Tahnke et al.ll2004D. W hile the stellar population diag- 
nostics of iKauffmann et al. (2003) are not immediately 
convertible into our simple colour indicators, their im- 
pressive and highly significant results point in exactly 
the same direction. 

While the results from all these redshift regimes are 
similar and point to a connection of nuclear activity and 
the presence of young stars, and that mass fractions of 
young stars are similar, we always find that a larger range 
of absolute masses is involved, showing as a range in UV 
luminosity. Here we find a variation of a factor of >40 
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for a given nuclear luminosity. 

Our host galaxy colours span a range that reaches the 
colours of Lyman Break Galaxies for a few very luminous 
hosts, while, as mentioned, the colours of the majority 
are somewhat redder than these. A comparison of op- 
tical/UV properties to the general population of high 
redshift galaxies would be very illuminating, but large 
statistical samples only become available in the near fu- 
ture, e.g. from the GOODS project. 

This persistent trend to find AGN to be associated 
with blue stellar colours is intriguing and suggests a close 
connection between enhanced star formation and nuclear 
activity. Additional support for such a connection comes 
from the detection of submm CO emission in a number of 
ext remely lumi nous high-redshift QSOs and radio galax- 
ies lOmontLet al. 2003), although the current sensitivity 
of submm telescopes is insufficient to perform this test 
for less luminous AGN at high z. 

While the fact that there is a relation can hardly be de- 
nied, its physical origin remains obscure. Is the enhance- 
ment of star formation a prerequisite for nuclear activity? 
Is it a simultaneously occurring phenomenon, caused by 
the same trigger? Or is it a consequence of the AGN? 
Galaxy merging and interaction are clearly two possible 
candidates to connect these two phenomena, but neither 
the only ones nor are the involved physics understood. 
Much additional data will be required, in particular those 
helping reliably to reconstruct the star formation history 
in high-redshift galaxies, before any firm conclusions can 
be drawn. 
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APPENDIX 

AGN IMAGES AND SURFACE BRIGHTNESS PLOTS OF AGN AND STARS 

Fig^ire 1 1 1 shows plots for each of the nine resolved object plus the co mp osite 'stacked' object. Two objects appear 
twice as they appear in overlapping areas of GEMS tiles. Figures IA12I and IA13I show a random selection of isolated 
stars used to show the zero case of a point sources without any host galaxy contribution, for comparison purposes. 
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Fig. All. — For each object the original images is shown and the uncorrected peak scaled PSF subtracted images, for both the F606W 
and F850LP filter. The image size is 2" X 2" and the images are shown in a linear grey scale, with overplotted logarithmic isophotes of 
0.5 dex spacing. The plots below show the radial surface brightness in the two filter, respectively. The curves show the data points with 
error bars from the original image (filled symbols), the PSF (dashed red line) and peak subtracted host profile (open symbols). The host 
galaxy is shown without any correction for oversubtraction as applied to the derived magnitudes. 
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Fig. All. — continued 



GEMS quasar host galaxies at 1.8 < z < 2.75 
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Fig. All. — continued 
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Fig. A12. — PSF peak subtraction applied to field stars: F 606W -band. Shown are surface brightness profiles for 12 randomly selected 
field stars in the F606W-band. Symbols and lines as in Figure IA 11 1 
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Fig. A13.— 



As Figure IA12I but for 12 different randomly selected stars in the F850LP-band. 



